Potential energy surface of the I'^A' Li2+Li doublet ground state 
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The lowest doublet electronic state for the lithium trimer [1^ A') is calculated for use in three-body 
scattering calculations using the valence electron FCI method with atomic cores represented using 
an effective core potential, ft is shown that an accurate description of core-valence correlation is 
necessary for accurate calculations of molecular bond lengths, frequencies and dissociation energies. 
Interpolation between 1^ A! ah initio surface data points in a sparse grid is done using the global 
interpolant moving least squares method with a smooth radial data cutoff function included in the 
fitting weights and bivariate polynomials as a basis set. The Jahn- Teller splitting of the 1^ E' surface 
into the 1^ A\ and 1^B2 states is investigated using a combination of both CASSCF and FCI levels 
of theory. Additionally, preliminary calculations of the l^A" surface are also presented using second 
order spin restricted open-shell M0ller-Plesset perturbation theory. 



Introduction 

With the success of ultracold molecular formation 
among the alkali metals over the past decade via 
photoassociation [1] and more recently with Fesh- 
bach resonances the dynamics of molecules in 
ultracold traps have become an important topic to 
many physicists. Alkali dimers have been formed in 
many different homonuclear and heteronuclear con- 
figurations, for both singlet and triplet ground elec- 
tronic states. Furthermore, recent experiments us- 
ing KRb O 0] and theoretical proposals for LiNa 
[5] have shown the possibility of efficiently forming 
ultracold ground vibrational state diatoms. While 
both ground and excited singlet and triplet states 
of alkali diatoms have been studied extensively both 
experimentally and theoretically, alkali trimers have 
generally been ignored by both theorists and experi- 
mentalists alike. In the last few years there has been 
an increase of interest in few-body physics with con- 
tinued success in the cooling and trapping of atoms. 

In these few-body systems, experimental interest 
in three body effects range from loss rate predic- 
tions to probing few-body quantum effects such as 
Efimov states (Esry and co-workers [B] and Grimm 
and co-workers [7] for example) due to the strong 
non-additive effects seen in alkali systems. To date, 
ah initio calculations for the sodium trimer have 
been done by several groups [5J [HI EH] as well as the 
potassium trimer [5J [TTJ [T^]. In the case of lithium, 
the quartet ground state I'^A' surface has been well 
studied [131 IHl 113 US] whereas the doublet system 
has been effectively ignored. The Jahn- Teller effect 
and conical intersections between the V^E' , and 
l^A' surfaces have been studied [17l [18], as well as 
the vibrational [19] , rovibrational [20] and hyperfine 
structures of the lithium trimer. However, to 
the best of our knowledge the entire ground or ex- 
cited state surfaces for the doublet lithium trimer 
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have not been completely studied. 

The structure of this paper is as follows, we first 
discuss the inclusion of core polarization potentials 
to accurately describe the effects of core valence 
correlation while using an effective core potential 
(ECP). The steps taken to optimize the basis set 
so as to provide accurate benchmark Li2 spectro- 
scopic values are then described in detail. Follow- 
ing this we show how accurate calculations of the 
potential energy surface (PES) were done with a 
low density of ab initio points through the use of 
the global interpolant moving least squares (IMLS) 
fitting procedure. Finally we describe the calcula- 
tions we have done on both the I'^A' and I'^A" sur- 
faces of the lithium trimer and future goals for the 
use of these surfaces. In this work, all calculations 
have been done using the MOLPRO 2008.1 quantum 
chemistry package [55] unless otherwise stated. 



Computational Details 

Core- Valence Correlation 

To accurately describe dissociation energies, equi- 
librium geometries and vibrational frequencies in 
alkali-metal clusters, it is necessary to account for 
the electronic core-valence (CV) correlation energy 
[23] . For all electron calculations this is possible for 
the lighter alkali atoms (Li through K) by using the 
explicitly correlated basis of Iron, Oren and Martin 
(lOM) 121]. This approach has recently been done 
by Cvitas et al [15 for the spin aligned I'^A' Lis 
surface using spin restricted coupled cluster calcu- 
lations with single, double and iterative triple ex- 
citations (RCCSD(T)). To account for CV correla- 
tion in heavier atoms, where all electron calculations 
are prohibitively expensive, it is necessary to use a 
core polarization potential (CPP). This is also nec- 
essary for valence electron full configuration interac- 
tion (FCI) calculations where an ECP (a physically 
equivalent representation to the frozen core approx- 
imation which has no CV correlation) is substituted 
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for the atomic electrons. This method has been used 
with great success for calculating both lithium dimer 
potential curves |25j and trimer potential surfaces of 
potassium |T2]. 

The theoretical description of the CPP is a 
straightforward addition to that of the ECP model 
of atomic cores. In the Born-Oppenheimer approx- 
imation the non-relativistic molecular Hamiltonian 
can be separated into kinetic and interaction opera- 
tors T and V respectively. Approximating the core 
of each nuclei with an / dependent pseudopotential 
and including the polarization effects at the nuclei 
gives for the interaction operator 

V = Y.^v^ + ^cpp) + E ;^ + (1) 

k j>i 

where 

^ nk ^ 

is the core pseudopotential, 

Pi'= = ^|fc?mi)(fc/m,| (3) 

m 

is the projection operator onto the subspace of an- 
gular momentum / on core k and Vcc is the core-core 
coulomb interaction. Here the polarization poten- 
tial, V^pp for a given core k, is expressed in terms of 
a static polarizability and external field at the nuclei 
position by 

Kpp = -^«feFfc.Ffc. (4) 

where the electric field at core k arising from the 
coulomb interaction with the electrons at v^i and 
other cores at R^j is 




The value of the static polarizability for Li+ is — 
0.1915ao [26] and the cutoS^ function C{rki) defined 
by Fuentealba et al [55] is given as 

C(rH) = (1 - e-^''^''--^ (6) 

with the cutoff parameter chosen to be 6k ~ 
0.831a(^^ [26]. This form of the cutoff parameter 
was first presented by Miiller ei al [57j and produces 
good agreement for ground and low excited states, 
however it does show diminished results for Rydberg 
states [57]. Our calculations using this core polar- 
ization potential to describe the core valence corre- 
lation energy were done using the MOLPRO 2008.1 
p2i implementation of the Fuentealba et al [26j CPP. 



TABLE I: Uncontracted basis set exponents for for the 
lithium Stevens, Basch and Krauss [28] pseudopotential 
basis used in this work. The s, p and d orbital exponents 
are each scaled to give an optimal dissociation energy for 
Li2 as discussed in the text. 



Orbital Type 


SBK 


LFK 


Scaled LFK 


s 


0.6177000 


0.6177000 


0.52504500 




0.1434000 


0.1434000 


0.12189000 




0.0504800 


0.0504800 


0.04290800 




0.0192300 


0.0192300 


0.01634550 


P 


0.6177000 


0.6177000 


0.00690155 




0.1434000 


0.1434000 


0.64858500 




0.0504800 


0.0504800 


0.15057000 




0.0192300 


0.0192300 


0.05300400 






0.0065729 


0.02019150 


d 




0.1346266 


0.13799227 






0.0546860 


0.05605315 






0.0180355 


0.01848639 






0.0076882 


0.00788041 



Basis Set 

Atomic and molecular polarizabilities are impor- 
tant factors for long range interatomic and molecular 
interactions. In ultracold systems this is a dominant 
contribution to the scattering length in addition to 
the location of the inner wall. To accurately describe 
both of these quantities in the lithium trimer, we re- 
quire that the basis set be flexible enough to describe 
atomic polarization at long range while accurately 
representing the inner wall. Very accurate polariz- 
ability calculations require an accurate description of 
electron correlation and a large basis set containing 
diffuse functions such as the Sadlej basis sets [31u32i. 
but this is computationally impractical to implement 
for large all electron systems. Accurate polarizabil- 
ity calculations can be achieved using an ECP with 
a small but well chosen basis set as demonstrated by 
Labello, Ferreira and Kurtz (LFK) [MUM]- Their 
basis set augments the Stevens, Basch and Krauss 
(SBK) ECP basis [28^, with an additional p function 
and four extra d functions optimized following the 
Sadlej [HI [S2 method. 

Using the SBK pseudopotential, we further opti- 
mized the uncontracted LFK basis (see Table [l] for 
the exponents) with and without the CPP for the 
Li2 X^T,'^ ground state. This was done using three 
scale factors p\{X = s,p,d) optimized at the FCI 
level to give the best calculated value of the dissoci- 
ation energy. As a benchmark, the CVQZ and CV5Z 
lOM basis p4] at the CCSD(T) level correlating all 
electrons (no frozen core) gives an error in the calcu- 
lated dissociation energy of 46.087 cm~^ and 27.902 
cm~^ respectively, while at the valence FCI level the 
unsealed LFK basis the error is 1.25 cm^^. Optimiz- 
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Method/Basis 


[A] 


ijJe [cm ^] 


Do [cm-^] 


[cm-'] 


Expt." 


2.673 


351.43 


8434.58 


8516.36'' 


Recommended'^ 


2.667 


352.98 


8340.12 


8516.43 


FCI/SBK scaled LFK 


2.693 


346.54 


8198.85 


8371.83 


FCI/SBK+CPP LFK 


2.663 


353.89 


8338.45 


8515.10 


FCI/SBK LFK 


2.687 


347.19 


8197.18 


8371.04 


RCCSD(T)/CVQZ'' 


2.676 


353.05 


8294.15 


8470.27 


RCCSD(T)/CV5Z" 


2.674 


353.08 


8311.70 


8488.46 



"Taken from Herzberg [29] 

''Extracted from the scattering RKR data |30| 

FCI/SBK+CPP with scaled LFK using Ps = 0.85, /3p = 
1.05 and /3a = 1.025 

CVnZ is the n-zeta core valence correlation consistent basis 
set from lOM [24]. 

TABLE II: Comparison of different spectroscopic constants of Li2 using benchmark basis sets with both explicit 
inclusion of core-valence correlation in RCCSD(T) calculations or through the empirical contribution through a core 
polarization potential with valence electron full configuration interaction. 



ing the P\ coefRcients to give the best dissociation 
energy, we obtain a difference of 0.077 cm~^ using 
the scaled exponents listed in Table |T] correspond- 
ing to the scale factors f3s — 0.85, f3p — 1.05 and 
(3d — 1.025. In Figure[l]and Table [TT| the results from 
the different basis functions and methods can be seen 
compared to the Rydberg-Klein-Rees (RKR) curve 
adjusted to reproduce the ultracold scattering re- 
sults ^35j. We chose to use this scaled, uncontracted 
LFK basis with the SBK ECP and Fuentealba CPP 
in calculating the potential energy surface of the 
lithium trimer at the FCI level of theory. 



Surface Representation 

Calculations of ab initio points at the FCI level 
are very computationally intensive, even for three 
electrons with the compact basis just described, so 
to accurately describe an entire potential energy sur- 
face with a low density of ab initio points we imple- 
mented the global IMLS fitting method 0^ . In this 
method the potential energy at an arbitrary point 
Z in the {x,y) plane is approximated by the use of 
a linearly independent basis bj(Z)(j = 1, . . . , n) and 
expansion coefficients aj{Z)(j ~ I,...,n) so that 
the interpolated energy is given by 



(7) 



The expansion coefficients a(Z) are found by mini- 
mizing the error function 



of the interpolated energy Vint and the ab initio en- 
ergy fi{ii){i — 1, • • . , Nd) at coordinates £i where Nd 
is the number of ab initio data points. 

Expressing the solution to the normal equations 
dE{Vint) / daj = in matrix form we obtain the lin- 
ear equation for the expansion coefficients 36J 



BW{Z)B'^Si{Z) = BW{Z){, 



(9) 



where W{Z) is the diagonal matrix of weights Wi{Z) 
and B is the matrix 



B = 



62(4) ■■■ 



hilNd) 



(10) 



The linear system in Eq. |9] is routinely ill condi- 
tioned and so is solved by singular value decomposi- 
tion. 

The weights Wi{Z) dictate the effective range at 
which a given ab initio point will contribute to the 
global fit and the effective contribution to the fit. 
We use Guo et al's [37 form of the weight func- 
tion, which introduces a cutoff function S{x) to the 
unnormalized weight function Ui(||Z — ^i||) so as to 
smoothly go to zero at a given cutoff radius Rcut- 
The cutoff function is given [37] by 

^(X) = l^^-^"^' ifO<X<l, 
^ ' [0 ifx>l, 

with m = 10 and the unnormalized weight function 
is 



(8) 



^{Z) = 



expH\Z-£,r/e] 
i\iZ-£A[/C)^ + e 



(12) 
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FIG. 1: Benchmark singlet Li2 potential energy curves using both the core-valence correlation consistent basis sets 
of Iron, Oren and Martin ;24- and core polarization potential Stevens, Basch and Krauss pseudopotential |28^ with 
an extended and scaled ECP basis set calculated at the RCCSD(T) and full configuration interaction level of theory 
respectively. The scattering RKR curve is the inner wall shifted version of Cote et al [30 such that the potential 
predicts the correct scattering lengths and Feschbach resonances. 




FIG. 2: Molecule frame Jacobi coordinates are used to 
describe the lithium trimer geometry. Within this coor- 
dinate system the potential energy surface is calculated 
assuming that the lithium dimer bond is adiabatically 
relaxed as the colliding lithium atom approaches. 



where the inverse coordinates Zj — 1/xi are used in 
this work. With the choice of coordinates used here 
there is a coordinate singularity in C2v symmetry. 
To avoid this all coordinates shifted by the 

same positive, arbitrary additive factor for the fit 
then transformed back upon completion. The scal- 
ing parameter ^ was chosen to be the average dis- 
tance between the interpolant point Z the ah initio 
points. With this definition the cutoff radius was 
defined in terms of as Rcut = 50.0 * C- This inter- 
polant method is used for the lithium trimer I'^A' 
PES to obtain a global fit using a low number of ah 
initio points as references. 



where e = 10 removes the singularity at ti. Then 
the normalized weight function is 

T.fS{\\Z-t,r/R,^t)v,{Z) 

where Rcut is to be determined as to give the best fit . 
Finally the basis functions are chosen to be bivariate 
polynomials of order n — Q such that 

h{Z) — \, Zi, Z2, Zi, Z2, Z1Z2T . ■ T Z^^ ^Z2,Z\Z2 ^, 

(14) 



Results and Discussion 

The ^^ A! surface of Lis was calculated at the full 
configuration interaction level using the scaled LFK 
basis set, the SBK pseudopotential [28 and core po- 
larization potential described above, with the three 
valence electrons included in the FCI calculation. At 
the FCI level, there are 410670 configurations of A' 
symmetry and 383292 configurations of A!' symme- 
try. All FCI calculations were performed with the 
Knowles-Handy determinant FCI program [SSJ 155] 
using the MOLPRO 2008.1 package 122,. The ge- 
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FIG. 3: Near equilibrium geometry potential energy surface for the Lis l^A' electronic state calculated using 
valence elctron full configuration interaction theory. Equilibrium is found to be at bond lengths of — 3.218A and 
Rc = 2.238A for C2v geometry configuration. 



ometry was chosen so to best describe the diatomic- 
atomic collision process. As such we used the molec- 
ular frame Jacobi coordinates for a homonuclear sys- 
tem where we define a vector Vf, along the diatomic 
inter-nuclear axis and another vector Rc from the 
diatomic center of mass to the colliding atom where 
the collision angle A is defined from the Coov axis 
(see Figure[2|. With this coordinate system the most 
efficient grid of ab initio points is an evenly spaced 
angular grid with the radial grid chosen to have the 
highest density of points at the minima. We cal- 
culated 26 with the collision angle varying between 
90° and 45° on the l^A' Lis surface by choosing 
A and Rc then optimizing the dimer bond length 
to give the lowest energy configuration. The l^A' 
state is found to have a triangular equilibrium ge- 
ometry on the l^Ai surface with rg = 3.218A and 
Rc = 2.238A for the Jacobi bond lengths with a 
dissotiation energy of 4979.42 cm~^. The interpo- 
lated surface near the equilibrium geometry is show 
in Figure [3j 

We have investigated the lowest states oi'^Ai, ^Bi 
and symmetry and the first excited state 
in point group C'2v using the complete active shell 
(CAS) method in conjunction with FCI. The same 
SBK and CPP representation of the core as in the 
FCI calculation of the l^A' surface was used here 
with an active space of 12 orbitals (5ai262462la2) 
for the energies of the l^Ai, l^Bi, l^B2 and 2'^Ai 
states as seen in Figure |4] Here Rc is fixed and rg is 
optimized at the CAS level with tight convergence. 
This is followed by a FCI calculation to obtain the 
energy at this geometry, with typical errors in the 
CAS geometry optimization compared to that of the 
FCI geometries on the order of a mA. A conical in- 
tersection between the l^Ai and l^i?2 surfaces is 
observed at Rc = 2.5lA which is the result of Jahn- 
Teller sphtting of the l^E' Djh surface lHj. 

Preliminary calculations of the Lis l^A" sur- 



CASSCK(3,12) Opumized HCI 




2.3 2.4 2,5 2.6 2.7 
RcLAJ 



FIG. 4: First four Lis doublet states in C2v symmetry 
near the l^Ai state equilibrium. A conical interesction 
is observed at Rc = 2.5lA between the l^Ai and 1^B2 
states with a further Lia minimum seen near Rc = 
2.25A for the l^B2 state. 



face have been carried out using second order 
spin restricted open-shell M0ller-Plesset perturba- 
tion (R0MP2) theory as implemented in Gaussian 



FIG. 5: The Lis l^A" potential energy surface calculated at the second order spin restricted open-shell M0ller- 
Plesset perturbation theory for collision angles 90° {C2v) to 45°. The dissociation energy is 14156. Scm"^ with 
Rc = 2.40A and re = 2.77A in the C2v symmetry. 



03 00]. Here the lOM CVTZ basis set H] was used 
for calculation size convenience, with a benchmark 
diatomic bond length error of 0.075A at R0MP2 
and O.OOSA at RCCSD(T). The analytic potential 
energy surface was interpolated using cubic splines 
with 256 ab initio data points and is shown in Fig- 
ure [s] This I'^A" state has a dissociation energy of 
14156.5 cm""'^ at the equilibrium triangular geome- 
try of Rc = 2.40A and = 2.11 k in C2v 



Conclusions 

The electronic ground state PA' surface of the 
lithium trimer has been calculated using valence 
electron FCI theory with the lithium cores repre- 
sented using the SBK pseudopotential [IS]. It was 
found to be necessary to systematically include core- 
valence correlation in the calculation for precise cal- 
culations of diatomic spectroscopic values. The ba- 
sis set chosen is a p and d function augmentation 
of the SBK basis set |2S] as given by Labello, Fer- 
reria, and Kuntz |33j with the s, p and d functions 
optimized with the inclusion of a core polarization 
potential to predict the correct diatomic dissocia- 
tion energy. With the recommended basis set of this 
work, the use of the core polarization potential to 
include core-valence effects lead to an improvement 
26.4 mA in the bond length, 6.44 cm^^ of the vi- 
brational frequency and 144.605 cm"^ dissociation 
energy of Li2 . 

To interpolate between sparse ah initio data 
points, the fuU interpolant moving least squares 



method was implemented using a scaled exponential 
weighting function with a smooth cutoff function as 
a constraint on the number of included data points. 
The surface was expanded using the inverse spa- 
tial coordinates with a 6 order bivariate polynomial. 
With this interpolation method the ^^ A! surface was 
calculated for collision angles 45° to 90° near the 
equilibrium Jacobi bond lengths of re = 3.218A and 
Rc = 2.238A. A conical intersection is found be- 
tween the l^Ai and 1^B2 states in C2v symmetry 
near the equilibrium geometry of the l^A' surface. 
Because of the location of this intersection it is nec- 
essary to account for its existence in both chemical 
and ultracold physics. Also a preliminary surface for 
the excited state I'^A" is presented at the R0MP2 
level of theory. It is the authors intention to continue 
to study the long range interactions of the lithium 
trimer on the ground l^A' surface and to investi- 
gate both elastic collisions and photoassociation of 
the lithium diatom-atom pair for the formation of 
ultracold trimers. 
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